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In this letter we report on a numerical investigation of the Aoki phase in the case of finite 
temperature which continues our former study at zero temperature. We have performed simulations 
with Wilson fermions ai P — 4.6 using lattices with temporal extension A^^ = 4. In contrast to 
the zero temperature case, the existence of an Aoki phase can be confirmed for a small range in k 
at /3 = 4.6, however, shifted slightly to lower k. Despite fine-tuning k we could not separate the 
thermal transition line from the Aoki phase. 
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I. INTRODUCTION 

The interest in the phase structure of Wilson fermions 
coupled to SU{3) gauge fields has revived recently, in 
particular due to the unexpected discovery of a first- 
order transition in twisted mass QCD This transition 
survives when the Wilson gauge action is replaced by a 
renormalization group improved gauge action 0. This 
has brought reports on a nontrivial phase structure 0,3 
back into the center of interest including the pioneering 
paper by Aoki Q. 

The Aoki phase, characterized by a non- vanishing con- 
densate (ipij^T^ip) and a broken pion mass triplet, was 
analyzed in the framework of chiral perturbation theory 
in Iq, . In [3 two scenarios were described: at a given 
gauge coupling (3 either the Aoki phase is realized in a 
certain k interval (separated by second order phase tran- 
sition lines from a phase with degenerate massive pions), 
or there are first order transitions in a certain interval 
of twisted masses (also ending at second order transition 
points) . 

In a recent paper Q we have found evidence that the 
Aoki phase at zero temperature is unlikely to extend to 
(3 > 4.6, whereas the mentioned new phase transition has 
been observed at /? = 5.2 The change between those 
two scenarios seems to happen between these two values 
of p. 

From the very beginning it was interesting to know 
where the Aoki phase lies at finite temperature. Numer- 
ical studies 0, have given some indications that the 
Aoki phase is completely embedded in the low tempera- 
ture phase as proposed by Aoki et al. |^. Their results 
for Nf = 2 and Nt- = 4 are shown in Fig. How- 
ever, the results left open the question, whether the Aoki 
phase will join the thermal transition line Kt{P, Nr) as 
it extends further in /? for a given thermal lattice exten- 
sion Nt- 

To be specific, the following observations were made by 
Aoki et al. for Nf = 2 dynamical flavors: For Nr = 4 
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FIG. 1: The phase diagram with Wilson fermions at finite 
temperature as presented in |^. We have added a point at 
/3 = 4.6 where we have found evidence for the Aoki phase. 



breaking condensate {'ipij^T'^il') was found in a k interval 
where m^± = 0. Surprisingly, the cusp of the Aoki phase 
moved slightly towards larger /3 when Nr is increased. 
In fact, going from N^ = 4 to Nr = 8 resulted only in a 
shift from /3 = 4.0 to f3 = 4.2. At larger f3 no evidence 
for the Aoki phase — at finite temperature — has been 
found in unquenched simulations so far. 

In a finite-temperature study by the MILC collabora- 
tion on a 12'^ X 6 lattice l3l| metastabilities of plaquette 
values at (/3,k) =(4.8,0.19), (5.02,0.18) and (5.22,0.17) 
were found. This observation could be compatible with 
the metastability at (5.20, 0.1715) observed in 0. 

A possible scenario could be that the tip of the cusp 
simply stops moving and turns into a first order line along 
which m^r ^ 0, separating phases with rriq < (at high 
k) from niq > (at low k) pTLIl2|. 

What happens to the first order phase transition line 
at weaker coupling? Is it universal and how does this 
affect the continuum limit? Does the transition extend to 
/3 ^ oo as a single first order line? Or does it split again 
giving room for a new Aoki-like phase enclosed? Can the 
Aoki phase extend towards /3 ^ oo being enclosed by 
two second order lines using improved actions? 
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FIG. 2: The real part of the Polyakov loop, its susceptibility 
XL and the order parameter (tpi'ysT'^tp) as a function of k at 
h = 0.005 and 0.001 fixed. The lower panels shows h = 0.003 
additionally. The data have been obtained from simulations 
on a 8^ X 4 and a 10^ x 4 lattice at /3 = 4.6. Throughout the 
same symbols are used. 
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FIG. 3: The (exponential) autocorrelation time r^xp (in terms 
of HMC trajectories) of the real part of the Polyakov loop is 
shown as a function of k for two different values of h. The 
lattice size is 10^ x 4. 



II. NUMERICAL RESULTS 

In this study we looked for the Aoki phase at finite 
temperature at /3 larger than in earlier studies. We 
performed numerical simulations at temporal lattice ex- 
tension N-r = 4 using dynamical unimproved Wilson 



fermions. The Wilson fermion matrix Mw was supple- 
mented by an explicit parity-flavor symmetry breaking 
term, i. e. the two-flavor fermion matrix was given by 



M (h) = Mw + h ij5T^ 



(1) 



We simulated at /3 = 4.6. k was fine-tuned in the in- 
terval [0.193,0.199] at fixed values of h. For each triple 
K, h) the order parameter {ipij^T^ip), the real part of 
Polyakov loop KeL and its susceptibility xl were mea- 
sured. These observables are shown in Fig. |21 From the 
Figure we see that ReL and Xh reveal a finite temperature 
transition around k = 0.19705, while the order parame- 
ter (ipij^T^ip) signals an existence of the Aoki phase in 
this K region. 

Note that we had found a vestigial region around 
K « 0.1984 at /3 = 4.6 in our previous zero-temperature 
study 0- Therefore, the Aoki phase seems to follow the 
finite-temperature transition line and is shifted to some- 
what lower K compared to the zero-temperature case. 
This shift, however, cannot be resolved in Fig.^ In any 
case, the Aoki phase extends inside a longer cusp than 
seen in P] . 

In the region of interest around k — 0.19705 we found 
rather long autocorrelations in the Monte Carlo time his- 
tories in contrast to other k values. The autocorrelation 
function has been estimated for the Polyakov loop and 
we found large values for the exponential Texp and the 
integrated autocorrelation time Tint (in terms of HMC 
trajectories). Examples for r^xp at /i = 0.005 and 0.001 
are shown in Fig. |21 

In order to study the Aoki phase, the limit 



lim lim {tpij^T tp) 



(2) 



has to be taken. Therefore, we made runs in the interval 
K e [0.1968, 0.19720] at several values of h e [0.001, 0.02] 
using spatial volumes 8'^ and lO'^. 

In the upper panels of Fig.^the results of those simu- 
lations are shown together with fits to the data from the 
largest lattice available at each h. For the fits we used 
the ansatz [a 



a{h) =A + Bh^ 



(3) 



where a = (ipij^r'^'ip). The value of the order param- 
eter at /i = is given by the fit parameter A. The fits 
are quite robust against the introduction of linear and 
quadratic corrections. Details of the statistics achieved 
are presented in Tab. ^ We did not list the information 
refering to the smaller lattice size 8^^ x 4 in all those cases, 
where 10^ x 4 data were available. The results of all our 
fits are quoted in Tab. ^ 

The lower panels of Fig. ^ show the same results as 
in the upper panels, however, in a different parametriza- 
tion (so called Fisher plots, see e.g. |^ ^3 This 
parametrization has the advantage that the curves bend 
upwards, if (ipij^T^ip) is non-zero in the limit h —^ 0. 
Hence, the finiteness of the intercept point can be read 
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FIG. 4: The upper panels show the order parameter {tpi^^r^ tp) as a function of h at five values of k. The lines are fits to the 
data as described in the text. The lower panels show Fisher plots of the order parameter and the corresponding fitting function 
at the same values of k as above. Throughout the same symbols are used. 
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TABLE I: Statistics used for the final analysis at selected k at /3 = 4.6. 



off better. In addition, it enables us to see directly if 
this intercept grows with the volume. At k = 0.19700 
we have computed the order parameter at lower values 
of h also for the smaller lattice size 8'^ x 4. As expected, 
a clear volume dependence is found, indicating that the 
existence of the Aoki phase becomes visible only for suf- 
ficient large volumes. The right most panels of Fig. 0] 
show data at (k = 0.19720) to illustrated a case where 
the order parameter {il^ij^T^ip) is found to vanish. 

III. DISCUSSION 

By measuring the order parameter (ij^ij^T^'ip) and ex- 
trapolating it to /i = we conclude that there is a parity- 
flavor broken phase at finite temperature [Nr = 4) at 



f3 = 4.6 in the interval 0.1968 < k < 0.19715. In the 
zero-temperature case we had found 131 that the Aoki 
phase seems to end in the vicinity of /3 = 4.6, k = 0.1984. 
In other words, for a finite temperature the endpoint of 
the Aoki phase is shifted towards larger (3 and also to 
lower K. 

This result differs from |^, ^3 where the endpoint 
of the Aoki phase for Nt = 4 was seen at /? = 4.0, 
K — 0.22. The finite-temperature transition line de- 
termined in Refs. 0, 0| crosses the interval f3 — 4.6, 
0.1968 < K < 0.19715 in which we see the Aoki phase. Al- 
though we have studied the respective k interval in rather 
small steps, we are not able to separate (the boundary 
of) the Aoki phase from the thermal phase transition 
(cf . Fig. . We plan to extend our investigation to larger 
Nr values in the near future. 
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TABLE II: The parameters of the ansatz a{h) = A + Bh'^ + Dh + Eh^ fitted to the data of {iiijsT^ip) at /? = 4.6 without 
(fits labeled 1 and 2), with linear, or with quadratic corrections (labeled 3 and 4). At each h the result from the largest lattice 
was used in the fit (for details see Table HJ. Fixed parameters are given by their value without indicating an error. In each case 
the first fit (bold numbers) was used in Fig. 2] 
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